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We consider a situation where the density and peculiar velocities in real space are 
linear, and we calculate the two-point correlation function in redshift space, incor- 
porating all non-linear effects which arise as a consequence of the map from real to 
redshift space. Our result is non-perturbative and it includes the effects of possible 
multi-streaming in redshift space. We find that the deviations from the predictions of 
, the linear redshift distortion analysis increase for the higher spherical harmonics of 

^ ' £^s- While the deviations are insignificant for the monopole the hexadecapole ^4 ex- 

, hibits large deviations from the linear predictions. For a COBE normalised F = 0.25, 

' /i = 0.5 CDM power spectrum our results for ^4 deviate from the linear predictions by 

a factor of two at the scales ~ 10/i~^Mpc. The deviations from the linear predictions 
' depend separately on /(fi) and h. This holds the possibility of removing the degener- 

acy that exists between these two parameters in the linear analysis of redshift surveys 
which yields only /3 = f{il)/b. 
, We also show that the commonly used phenomenological model where the non- 

C^h' linear redshift two-point correlation function is calculated by convolving the linear 

Q , redshift correlation function with an isotropic pair velocity distribution function is a 

^ ' limiting case of our result. 

^ I Key words: Keywords: Cosmology: theory — cosmology: observation — dark matter 

— galaxies: distances and redshifts — galaxies: clustering — large-scale structure of 
universe 



1 INTRODUCTION 

It has been long recognized that galaxy peculiar motions introduce distortions in the clustering pattern observed in redshift 
surveys. On small scales, the random motions of galaxies in virialized clusters causes structures to appear elongated along the 
line of sight, often referred to as the Fingers-of-God. On large scales, the coherent inflow of galaxies into over dense structures 
causes these to appear flattened along the line of sight (Peebles 1980, §76, and references therein). 

In a seminal paper Kaiser (1987) showed that in the linear regime the galaxy peculiar velocities are expected to introduce 
a quadrupole anisotropy in -Ps(k) the power-spectrum in galaxy redshift surveys, the relation between Ps(k) and its real space 
counterpart Pr{k) being 

P.(k) = (l + f^EJV')'^'r(fc) (1) 

where ^ is the cosine of the angle between k and the line of sight. He proposed that this effect could be used to measure the 
cosmological density parameter Q. This possibility has led to a lot of work on Linear Redshift Distortions and the reader is 
referred to Hamilton (1998) for a comprehensive review. Although in principle Kaiser's original proposal is quite simple, there 
are a large number of problems which arise in the actual analysis of redshift distortions. 

One of the main problems arises from the possibility that galaxies may be biased tracers of the underlying mass distribution 
which determines the peculiar velocities. The fact that galaxies of different types cluster differently (e.g., Dressier 1980; Lahav, 
Nemiroff & Piran 1990, Santiago & Strauss 1992; Loveday et al. 1995; Hermit et al. 1996; Guzzo et al. 1997) implies that not 
all of them can be exact indicators of the mass distribution. The presence of bias is a consequence of the complex process 
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of galaxy formation and this suggests a scale dependent, stochastic and nonlinear bias (e.g. Cen & Ostriker 1992, Mann, 
Peacock & Heavens 1998, Ue-Li Pen 1998, Dekel & Lahav 1999). Various studies (e.g. Taruya et al. 2001, Benson et al. 
2000, Kauffmann, Nusser & Steinmetz 1997) indicate that deterministic linear biasing is a reasonably good assumption at 
large scales where the density fluctuations are linear. In this situation biasing is described using only one scale independent 
parameter b which relates the fluctuation in the galaxy density and the mass density through a linear relation 5g — bSm- It 
then follows that the analysis of linear redshift distortions allows the determination of the quantity /3 = O!^'^ /b, and not Q and 
b individually. The analysis of different redshift surveys has given values of (3 which vary considerably in the range 0.2 — 1.1 
(Table. 1 of Hamilton 1998) . A part of the large spread in values is related to uncertainties in the nature and extent of bias 
in the different galaxy samples. The parameter /3 can also be measured by comparing the observed peculiar velocities with 
the observed galaxy distribution (eg. Strauss and Willick, 1995). The spread in the values of (3 determined using this method 
(Table 3. of Strauss and Willick, 1995) is comparable to that in the values determined using redshift distortion. Another 
method of measuring /3 is based on the observed cluster abundance (Wu, 2000). Melott et al. (1998) have proposed that the 
"Bulls-eye effect" which arises due to redshift distortion offers a different way of probing fl. 

Kaiser's original work and much of the subsequent work is based on the plane-parallel-approximation (PPA) which is 
valid if the angles subtended by the pairs of galaxies in the analysis are small. The radial nature of the redshift distortions has 
to be taken into account when analyzing wide angle surveys. Different strategies for incorporating this effect have been studied 
by Fisher, Scharf & Lahav (1994), Heavens & Taylor (1995), Hamilton & Culhane (1996), Zaroubi & Hoffman (1996), Szalay, 
Matsubara & Landy (1998) and Bharadwaj (1999). The radial gradients of the galaxy selection function and the motion of 
the observer introduce new effects when the radial nature of the redshift distortions are taken into account. These two effects 
do not contribute in the plane parallel approximation. 

The linear analysis of redshift distortions, applicable on large scales where the real space density is linear (a; < 1), 
incorporates the distortion due to peculiar velocities only to linear order when transforming from real to redshift space. N- 
body simulations suggest that non-linear effect are significant in redshift space even at considerably large scales where the real 
space density field is linear (e.g. Suto & Suginohara 1991; Fisher et al. 1993; Graman, Cen & Bahcall, 1993; Brainerd et al., 
1996; Bromley, Warren & Zurek, 1997). This is also seen in the two-point correlation function measured in redshift surveys. 
In addition to the flattening predicted by the linear analysis, the two-point correlation function remains elongated along the 
line of sight at scales as large as 20^~^Mpc (e.g. Figure 2. Peacock et al. 2001). This indicates that there is a regime where 
the density and peculiar velocity fields in real space are adequately described by linear theory, and it is the mapping from 
real to redshift space which is non-linear. This is often referred to as the translinear regime. It is important to incorporate 
the non-linear effects in the analysis of redshift distortions in this regime and there have been two different approaches to this 
problem. 

The first approach (e.g. Fisher et al. 1994, Peacock & Dodds 1994, Ballinger, Peacock & Heavens 1996) ) is phenomeno- 
logical and the effect of the non-linearity is incorporated by convolving the linear redshift two-point correlation function with 
the line of sight component of a random isotropic pairwise velocity distribution function (eq. ^9|). Equivalently, these effects 
are introduced into the power spectrum through a function f{k,fj) which multiplies the linear redshift power spectrum. The 
resulting redshift power spectrum is given by 

p.(k) = {l + p^l^fPr{k)f(k,^JL) . (2) 

Different forms, including a Gaussian and an exponential, have been used for the distribution function and this approach is 
found to match the results of N-body simulations well (Hatton & Cole 1998). 

The other approach is based on the Zel'dovich approximation. Taylor & Hamilton (1996), Fisher & Nusser (1996) and 
Hatton and Cole (1998) have used the Zel'dovich approximation to analytically study the behaviour of the redshift-space 
power spectrum in the translinear regime. They find that the results from Zel'dovich approximation are in agreement with 
N-body simulations in predicting the shape of the quadrupole to monopole ratio for the redshift power spectrum. They also 
demonstrate that there can be departures from the linear predictions even at relatively large scales. The first two groups 
of authors also conclude that the agreement between the Zel'dovich approximation and N-body simulations suggests that 
it is the coherent infall into over dense structures and not random motions in virialized clusters which is responsible for 
departures from the linear behaviour in the translinear regime. Recently Hui, Kofman & Shandarin (2000) have used the 
Zel'dovich approximation to study the probability distribution function of density in redshift space. They find that the 
peculiar velocities may significantly increase the multi-streaming in redshift space even when this is insignificant in real space. 
This suggests that non-perturbative effects may be important in the mapping from real to redshift space even when the real 
space density is relatively linear. 

It is the aim of this paper to analytically study the redshift two-point correlation function in the translinear regime 
focusing on the non-linear effects which arise from the mapping from real to redshift space. It is hoped that such a study 
will help elucidate the relation between the redshift correlation function and the different real space quantities on which 
it depends. We consider a situation where the density fluctuations and peculiar velocities are linear in real space, and we 
calculate the two-point correlation function in redshift space taking into account all non-linear effects of the redshift 
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distortion. The calculation closely follows the method used to calculate the real space two-point correlation function in the 
Zel'dovich approximation (Bharadwaj 1996b). Our result is non-perturbative and hence it incorporates the effects of multi- 
streaming in redshift space. Our calculation assumes the plane-parallel approximation and hence the effect of gradients in the 
selection function and the motion of the observer have been ignored. The calculation is presented in section §2. 

In §3 we compare our results with the predictions of linear redshift distortion and investigate the nature and extent of the 
deviations from the linear results. The effects of linear redshift distortion are quantified by a single parameter /3 = f{Q)/b, and 
it not possible to determine fl and b separately. We investigate whether this degeneracy is broken and if it becomes possible 
to determine fl and b if the non- linear redshift distortions are taken into account. 

In §4 we present the summary and discussion. 



2 CALCULATING 

The position of a galaxy in redshift space si differs from its actual position ri, the relation between the two being 

si - Uini = ri (3) 

where ?7i = vi • ni is the line-of-sight component of the peculiar velocity of the galaxy and units have been chosen so that the 
Hubble parameter H = 1. We study how the galaxy two-point correlation function in phase space p2r(i'i, r2, vi, V2) is related 
to its redshift space counterpart p2s(si, S2, vi, V2). 

The mapping from real space to redshift space preserves the number of galaxies 

P2s{si,S2,-vi,-V2)d^sid^S2d^vid^V2 = p2r(ri,r2, vi, V2)d^rid^r2d^wid^i'2 (4) 

For pairs of galaxies which are at a large distance from the observer and subtends a small angle in the sky, the Jacobian of 
the transformation from s to r (eq. ^) may be neglected and we get 

P2s(si,S2,Vl,V2) = P2r(si - (7in, S2 - f/2n, Vi,V2) (5) 

in the "plane parallel approximation". Here the unit vector n refers to the common line-of-sight to the pair of galaxies, and 
the line-of-sight components of the peculiar velocity of the two galaxies are taken to be parallel. 

The phase space two-point correlation function in real space is homogeneous which allows us to write equation (^) as 

P2s(si,S2,Vl,V2) = P2r(s - f/n, Vl,V2) (6) 

where s = S2 — si is the vector joining the pair of galaxies and U = n - (v2 — vi) is the line-of-sight component of the difference 
in peculiar velocities of the two galaxies. The vector s is decomposed into two parts 

s = sx + S||n, (7) 

S|| the component along the line-of-sight and sx the part perpendicular to the line-of-sight. 
A Taylor expansion of equation ^ in powers of U gives us 

P24Sl,S2,Vi,V2) = 2^ ^1 g^P2r(s,Vi,V2). (8) 

We calculate ^s{s±, sy) the galaxy two-point correlation function in redshift space by integrating out the velocity information 
in equation 

r 

^sfl + Cs(s±, S||)] = / d^Vld^V2p2s{s,-Vl,V2) 

dVd^«2 5]^^^P2.(s,Vi,V2). (9) 

where Ug is the mean number density of galaxies and s± —\ s± |. Equation (^) relates the galaxy two-point correlation in 
redshift space to a sum of velocity moments of different orders in real space. The velocity moments in real space can also be 
expressed as 

d\ld\2U"'p2r{s,Vl,V2) = ^^{[/'"(l + + fe-)) (10) 

where {) denotes ensemble average, and 5i and 82 refer to the perturbation in the galaxy number density at the position si 
and S2 in real space. Using this and separating the odd and even powers of U in equation we obtain 

l + ?.(«x,S||) = E(2^(9^j {t/^'(l + M(l + M) 
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- E(^(^) (^--(i + M(i + ^..)). (11) 

q — 

Our analysis is restricted to a situation where the peculiar velocities and density perturbations are linear in real space. 
It is also assumed that these quantities are a Gaussian random field. We next discuss some of the statistical properties of U 
and S in real space. 

The ensemble average of an odd number of U's and S's is zero. For example 

(U) ^ (Sir) = (USlrSir) ^ . (12) 

The quantities in terms of which we express all the velocity moments encountered in equation (^) are the real space 
galaxy two-point correlation 

C.(sx,S||) = {5i52), (13) 
the line-of-sight component of the pair velocity 

Vp{s^, S||) = {U{5i + 52)) = n ■ {(V2 - vi)(5i + 82)) , (14) 
and the dispersion of the line-of-sight component of the pair velocity 

cr|>(sx,S||) = (U^) = ninj((v2 - vi)i(v2 - vi)j) . (15) 
We also have the relation 

Vp{s^,S^\)^2{U5lr)^2{U52r). (16) 

which follows from the statistical isotropy and homogeneity of the peculiar velocities and perturbations. 

All the higher velocity moments can be expressed in terms of ^s, Vp and ct|>. In addition, these three quantities are not 
independent and they can be related using linear perturbation theory (§4). 

We first consider the term {U^'') which can be expressed as 

{U''')=^j^[aUs.,s,)V. (17) 
We next consider the even velocity moments in equation (^t|). These can be expressed as (for g > 0) 

{U^''{l + 5lr){l + S2r)) ^ (18) 
= ((1 + 5l.)(l + S2r)) + {2q){2q - 1){U^^~^) (f/5l,.) (USir) 

= ■^kp(5±,g||)]^[l + er(s±,S||)]+ ^^_^^^_^^, [ap(gx,g||)]^ , 

and for the odd moments we have 

l^U^1+^{l + 5^,){l+62r)) = {2q+l){U^^){U{S^r + 52r)) 

= ^%^kp(sx,s||)]«l^p(sx,s||). (19) 
Using these for the odd and even terms in equation (|l^) we get 

l + C.(sx,S|i) = ^ _([a|,(sx,S||)]''[l + ?.(sx,S||)] (20) 

9 = ■ ''11 

~ a7j^{[o-p(sx,S||)]'"l/p(sx,S||)}+ ^{kp(sx,S||)]« } 

Each of the three terms in the above equation is next expressed using a Dirac Delta function. 



'''^^ " ^' ^ " j ^""^^'^^ " ^' ) ■ ^] ■ (21) 

Below we show this procedure for the third term in equation (po|). 

A[(Jp{s)Y ^} = d s S (s - s ){[ap{s )Y ^} 



Q^2,+ 2 LL J V /J 4 J ^^2,+2 

d'^'^^ / (C^exp[i(s-s ).k]hfcfa|,(s )]^ (22) 
where fcy = k ■ n is the line of sight component of k. 
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Carrying out a similar procedure for all the terms in (EOl) we obtain 



l+Ca(s±,S||) = 



Cr(Si,S||)+ 1- 



2 as II 



/d^k .,, \ , 1 1 r , 2 2 / ' ' 

^^exp[j(s-s )-k]2^— ap(sx,S|| 

9=0 

The sum over q gives a Gaussian in fcii which can be integrated 



(27r)3 



exp[i(s — s ) ■ k] exp 



= 5^(sx - Sx)G(s|| - S||,o-p(sx,S||)) , 



where we use 



(23) 



(24) 



(25) 



to represent a normalised Gaussian distribution. 

Using this in equation ( p^ we can integrate over d^sx to get 

2n 



1 +Cs(sX:S||) ^ ds 



Cr(sx,S||)+ 1- 



Vp(sx,S||) a 



G(S|| - S||,ap(sx,S||)) 



(26) 



2 as II 

Replacing d/ds^^ with the derivatives of the Gaussian distribution and making a change of variable S|| — Sy gives us 

1+C4s±,si|) = / ds|| G(s||,(jp(sx,S|| + sj,)) X (27) 



Cr(sx,S|| +S||) + 1 - 



S||Vp(sx, S|| + S||) \ V|(sx, S|| + S||) 



2crp(sx, s II + s!|) / 4crp(sx, s II + s!. 



This expresses the galaxy two-point correlation function in redshift space as a function of there quantities in real space, 
namely the galaxy two-point correlation function ^r, the line-of-sight component of the pair velocity Vp and the dispersion of 
line-of-sight component of the pair velocity ap. Equation (|^ is non-perturbative and it incorporates all the effects of redshift 
space distortions. This includes multi-streaming in redshift space and all other non- linear effects which may arise due to the 
mapping from real to redshift space. 



3 COMPARISON WITH LINEAR RESULTS. 



Kaiser (1987) and Hamilton (1992) have analyzed the effect of linear redshift distortions on the two-point correlation function. 
We use to refer to this and in this section we compare ^f" with which incorporates non-linear redshift distortions. 
The expression for is obtained by keeping only the linear terms in equation (EB) which gives 



d 



1 



(SX,S||) = ?r(sx,S||) - — yp(sx,S||) + 7; 7^2 P ' ^ \\) 



ds 



2 dsl 



(28) 



In the linear regime, under the assumption of linear bias, the real space quantities ^r, Vp and ap can be expressed in terms 
of ^, the two-point correlation function of the underlying dark matter distribution (details in Appendix A), as 

C.(sx,S||) = b'C(s) (29) 

where b is the linear bias parameter and s = ^ /s-^ + . We also have 



Vp(sx,S|| 
and 

crp(sx, S|| 



(30) 



= fm 



15 



Us) 



(31) 
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where /(fi) ~ fi'^'^ is the dimensionless growth rate for hnear perturbation (Peebles 1980), is the one-dimensional peculiar 
velocity dispersion and 

Us) = ^ r ay)y''dy. (32) 

* Jo 

We use equations (p9|),(^o|) and ( |3l| ) to express both (eq. and (eq. ^) in terms of only three inputs namely f{Q,),b 
and ^. This gives an operator equation (eq. 4 of Hamilton 1992) 



Cf(sx,S||) 



2 -12 



as) (33) 



for the linear redshift two-point correlation in terms of ^ the actual two-point correlation of the underlying dark matter 
distribution. 

An alternative way to parametrize the anisotropy of the redshift two-point correlation function is to use s and = sy/s 
(the cosine of the angle between s and the line of sight n) instead of {s±, S||). The angular dependence of ^s{s, /i) can be very 
conveniently expressed as a sum of spherical harmonics (Hamilton 1992) as 

(s) = Co (s)Po(m) + ^2{s)P2{fi) + ?f (s)P4(m) (34) 

where Piifj) are the Legendre Polynomials and i,Q , (,2 ans ^4 are the monopole, quadrupole and the hexadecapole components 
of the linear redshift two-point correlation function. In the linear analysis these 

Co^(s)=6^(l + |/3+i/32)C(s) (35) 
^Hs)=b\^(3+^f3'ms)~Us)) (36) 

(«) = m + f 6(s) - lUs)]; (37) 

are the only non-zero harmonics, and here l3 = f/b. The two-point correlation function of the underlying dark matter 
distribution is largely undetermined and equation (^Hj) cannot be used to determine both b and /3. The ratios of £,0 ,£,2 a-nd ^4 
can be used to only determine /3. It is a limitation of the linear analysis that b and /(fi) cannot be determined individually. 
The angular dependence of ^s(s,/i) can also be expanded in terms of spherical harmonics 

00 

1=0 

and here, as in the linear analysis, the odd terms are all zero. However in this case there are no simple analytic expressions 
for the spherical harmonics and these have to numerically evaluated. 

We have used a COBE normalised (Bunn and White, 1996) CDM power spectrum (Efstathiou, Bond and White, 1992) 
with shape parameter F = 0.25 and h = 0.5 to calculate the real space two-point correlation function for the underlying 
dark matter distribution. This model predicts CT8 = 0.55 and we expect non-linear effects to be small in real space at scales 
> 10/i~^Mpc. Below we compare the spherical harmonics of £s{s,fi) and £,3(3,11) for this model for different values of / and 
b. 



We first consider the case / = 1, 6 = 1 for which the various spherical harmonics are shown in figure Bl. We find that 
the deviations from the linear predictions are more pronounced for the higher spherical harmonics. While the behaviour of 
the monopole £q shows practically no deviations from ^q', ^4 exhibits significant deviations from s £4 . Another feature is that 
the non-linear redshift distortions produce a non-zero £e whereas the linear analysis predicts a value of zero for this and all 
higher harmonics. We have not considered any of the higher harmonics here. 

We next analyze the behaviour of ^2 and ^4 in some more detail using the quantity Ri{s) — £i{s) /£i' (s) whose deviation 
from the value Ri{s) — 1 indicates departures from the linear predictions. The behaviour of R2 and R4 are shown in figures 



B2 and 



B5 . We see that the behaviour of R2 and R4 looks similar, the difference being that the magnitude of the deviations 
from the linear predictions is much larger for the latter. A feature common to R2 and R4 is that the deviations from the linear 
predictions increase monotonically as we go to smaller length-scales for low f{fl), , whereas for larger f(fl) the deviations 
saturate and even fall at smaller scales. A possible explanation is that at small scales the anisotropics in are erased by 
multi-streaming in redshift space. The peculiar velocities increase with / (51) and we can expect the effects of multi-streaming 
to also increase with /(f2). We find that ^4 is around a factor of two larger than £4 on scales 10 — 20/i~^Mpc where the 
non-linear effects are expected to be small in real space. The value of R4 increases at scales smaller than 10/i~^Mpc but real 
space non-linear effects are also expected to become significant at these scales and our results are not expected to give the 
true picture here. Both R2 and R4 approach one at large scales. 

The b and /(fl) dependence of the non-linear redshift space distortions cannot be expressed through just one parameter 
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p. We consider the behaviour of the ratio €,4{s) / £,o{s) (figure B4) to study this effect. Linear redshift distortions (eqs. ^ and 
|3^ ) predict this to depend on just (3. We see that the our predictions for f4(s)/fo(s) change considerably if we vary /(fl) 
and b keeping (3 fixed. This holds the possibility of using the redshift two-point correlation function to determine f(yi) and b 
separately. 



4 DISCUSSION AND SUMMARY. 

We have calculated the galaxy two-point correlation function in redshift space for a situation where the density fluctuations 
and peculiar velocities are linear in real space. Our calculation incorporates all non-linear effects which arise due to the mapping 
from real to redshift space. Our result (equation is non-perturbative and it includes the effects of multi-streaming in redshift 
space. 

The position of a galaxy in redshift space is a combination of its actual position and the line of sight component of its 
peculiar velocity, and the redshift two-point correlation function will have contributions from three effects 

1. Correlations between the actual positions of the galaxies. This is quantified by the real space two-point correlation £,r{s). 

2. Correlations between the peculiar velocities and the actual positions of galaxies. This is quantified by the mean pair 



velocity whose component along the light of sight is Vp(s). In the linear regime this is negative (eq A6) because of coherent 
fiows out of under dense regions and into over dense regions. 

3. Correlations between the peculiar velocities of galaxies. This is quantified by the pair velocity dispersion whose component 
along the light of sight is crp(s) = 2[(t^ — (n- v(si)n- v(s2))] where s = S2 — si. Here is the one dimensional peculiar velocity 



dispersion of the random motion of individual galaxies and (n • v(si)n- v(s2)) — — /^■^-2-(V^) ^£,{s) (eq. A8) is the correlation 

d3|| I Y 

in the line of sight component of peculiar velocities due to coherent flows. At large separations the effect of the coherent 
fiows is much smaller than the contribution from random motions and (t|>(s) ~ 2(t^. The velocity- velocity correlations due 
to coherent fiows increases at smaller separations. This causes crp(s) to be anisotropic and its value is less than 2a^ . In the 
linear analysis cr|>(s) decreases monotonically as s is reduced and (j|>(0) = 0. 

The linear redshift two-point correlation (eq. incorporates these three effects and S,s{s) at a separation s in redshift 
space is expressed in terms of ^r(s),yp(s) and (Tp(s) at the same separation in real space. The line of sight component of the 
relative peculiar velocity U between a pair of galaxies causes their separation in redshift space to be different from the actual 
separation, and in principle the real space quantities should be evaluated at a different separation r where ry = s\\ — U and 
rx = sx. A possible way to incorporate this effect is to use (eq. for ^f'(s) and evaluate all the real space quantities at r, 
assuming that the different values of U being given by a probability distribution fiU). This gives 



= J dUfiU)eis^,s^i^U). (39) 

The distribution of U is characterized by o-p(sx,S||) which in general is anisotropic and scale dependent. As a simplifying 
assumption the distribution of U is usually taken to be isotropic and the function f{U) is chosen to be either a Gaussian or 
an exponential with only one constant parameter 2a^ which is the value of (Tp(s) at large separations. This phenomenological 
model for the non-linear redshift two-point correlation has been found to match the results of N-body simulations. Two effects 
which are not included in this model are (1.) the anisotropy and spatial dependence of crp(s) arising from the coherent fiows, 
and (2.) the fact that the distribution of U is correlated with the density fluctuations. 

Our exact calculation takes into account all of the effects discussed above. We find that the redshift two-point correlation 
function (eq. can be expressed as the real space two-point correlation function combined with ratios of the pair velocity 
to the pair velocity dispersion, convolved along the line of sight with a Gaussian one dimensional pair velocity distribution. 
To get a better understanding of the various terms in (eq. we first consider a limiting situation where the separation s is 
taken to be large so that (n ■ v(si)n ■ v(s2)) i.e. the contribution to crp(s) from the random motions is much larger than 

the contribution from coherent fiows. In this situation we can ignore the spatial variation of crp(s) (except for the contribution 
of the derivative of (Tp(s) to ^^(s)). The calculation of in this limit is presented in Appendix B. We find that in this 



limit the result of our calculation (eq. 84) is identical to the phenomenological model (eq. ^9|) with a Gaussian pair velocity 
distribution. This can be interpreted in terms of a simple picture where all the non-linear effects of the redshift distortions on 
arise from the random motions of galaxies. The random rearrangement of the galaxy positions over the length scale a along 
the line of sight in redshift surveys modulates the linear redshift correlation function (s). This diffusion process acting on 
5f'(s) along the line of sight introduces an elongation in ^s. 

It is interesting to note that a similar interpretation is also possible for the non-linear effects on the two-point correlation 
function in the Zel'dovich approximation (Bharadwaj 1996b) and in the perturbative treatment of gravitational instability 
(Bharadwaj 1996a). In both situations the scaling properties of the non- linear effects is determined by the pair velocity 
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dispersion (Bharadwaj 1997) and at large scales the non-linear two-point correlation function can be interpreted in terms of 
the linear two-point correlation function modified by small scale random motions. 

We now come back to the interpretation of equation (^^. One of the main features is that it incorporates all the effects 
of the anisotropy and scale dependence of the pair velocity dispersion. As noted earlier, these two features arise because of 
the coherent flows. Whereas the convolution in equation ( |39[ ) accounts only for the effects of random motions, the convolution 
with the distribution function in equation also includes the effect of the coherent flows. The coherent flows reduce the 
width of the distribution function, causing a reduction in the elongation in . The effect of the coherent flows is also present 
in the terms involving ratios of Vp and ap, the contribution of these terms increasing with the effect of coherent flows. All 
the terms involving Vp arise due to the non-linear redshift distortions. 

In the linear analysis the redshift distortions produced by the coherent flows give rise to a flattening of the two-point 
correlation function along the line of sight. This effect is best captured by the quadrupole moment ^2 which is predicted to 
be negative. The linear analysis also predicts a positive hexadecapole moment arising from the term involving the second 
derivative of (Tp(s) in equation ( p8| ) for ^f^(s). 

Comparing our result with the predictions of the linear analysis of redshift distortions we see that the non-linear redshift 
space effects are more pronounced on the higher spherical moments of the redshift space two-point correlation function. 
While the monopole ^0 shows no significant deviations from the linear predictions, the deviations are around 10% on scales 
of 10 — 20/i~^Mpc for ^2. The moment ^4 shows very large deviations from the linear predictions, the non-linear predictions 
differing by a factor of ~ 2 on scales of 10 — 20/i~^Mpc. 

A point which should be noted is that our analysis ignores all non-linear effect in the real space density fiuctuations and 
peculiar velocities, and focuses only on the non-linear redshift distortions. While this assumption is expected to be valid for 
and Vp at scales greater than 10/i~^Mpc, we do not expect this to hold for ap which has a constant isotropic part 2a^ arising 
from the random motions. This has a large contribution from small scales which are non-linear and we do not expect the linear 
results used here to give a realistic estimate of a%. We also expect this effect to be present in our estimates on the effects 
on non- linear redshift distortions on the different angular moments of ^s. A possible solution is to use N-body simulations to 
determine the constant, isotropic part of ap and to use linear theory to calculate the contribution of the coherent flows. This 
possibility has not been explored here, and we plan to continue future work in this direction. 

Another possibility which might arise but has not been taken into account is non-linear bias. A quadratic bias relation 
gives rise to two bias parameters 61 and &2, and the relation between the fluctuations in the galaxy number density and the 
matter density is 5g — biSm + f 62<5m- In this situation 5g will no longer be a Gaussian random fleld and this will give rise to 
a new series of terms in equation (^o|) which could be summed up following the same procedure. 

The effects of linear redshift distortion can be parametrized by a single parameter /3 — f{Q,)/b and the analysis of redshift 
surveys based on this do not yield /(H) and b separately. We flnd that ^4, the hexadecapole moment of ^s, depends on /(H) and 
6 separately when the non-linear effects of redshift distortions are taken into account. This holds the possibility of separately 
determining /(SI) and b from redshift surveys. Although this effect increases at small scales, the range 10— 20/i~^Mpc is possibly 
best suited for putting it to use. Non-linear effects in real space and non-linear biasing will become important at smaller scales. 
At these scales the hexadecapole moment is expected to be 10% of the quadrupole moment. Further investigations are needed 
to determine the range of validity of the results presented here before they can be applied to the analysis of redshift surveys. 
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APPENDIX A: 

In this appendix we discuss, in some detail, the relation between ^r{r), Vp(r), a%{r) and ^(r). In the linear regime S{ri), the 
fluctuation in the dark matter density, can be expressed in terms of a potential i/'(ri) as 

5(ri) = V^^(ri), (Al) 
and the peculiar velocity v(ri) can be expressed as 

v(ri) = -/(S7)VV'(ri), (A2) 
In addition, assuming a linear bias relation, (5g(ri) the fluctuation in the galaxy number density can be expressed as 

<5s(ri) = 6V'^(ri). (A3) 
Using these we can express ^(r) as 
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CM = {5(ri)5(r2)) = {VV(ri)V^^A(r2)) 



(A4) 



where r = r2 — ri and </>(r-) = {tp{ri)tp{r2)) ■ Equation (A4) can be inverted to express (t>{r) and its derivatives in terms of 
We use these to calculate Vp{r), which gives us 



Vpir) = {n-(v(r2)-v(ri))(<5s(r2)+5(ri))) 

= (-&/){n ■ (V^/-(r2) - V7A(ri))(VV(r2) + VV(ri))) 
= -2/6(n- V)VV(r) 



(A5) 



This can be further simplified using equation (^4) to obtain 
Vp{r) = -2fb^{\7Y'ar) 



6(0 



(A6) 



We next consider cr|>(r) which can be expressed as 

'^p(r) = (n- (v(r2) - v(ri)) n- (v(r2) - v(ri))) 
= -2/ < [n- ^(rOlin- Vi/.(r2)])+a" 



(A7) 



where (t^ (n- v)'^ > is the one dimensional peculiar velocity dispersion. The expression for crp(r) can be further simplified 
using (A4) to obtain 

crp(r) = 



= f 



(A8) 



APPENDIX B: 



In this appendix we consider a limiting situation where the separation s at which the redshift two-point correlation is being 
calculated is very large. We also assume that the real space two-point correlation function is a power law ^r{s) oc x~'' with 
7 > 2 at large separations. The power law requirement is not crucial and the arguments presented here are valid provided the 
integral Cr(s)sds converges for x oo. 



With the power law assumption we have Vp{s) oc ^ (eq. A6), o-|>(s 



(eq. A8) and (Jp(s) ~ 2a for large 



values of s. The point to note is that at large separations the correlation between the galaxy peculiar velocities is very small 
and the pair velocity dispersion is nearly isotropic and has a constant value 2a^. 

We next shift our attention to equation (^) which expresses as a series of different terms involving £,r,Vp and a%. In 
the linear regime each of these quantities is characterized by a small number e ~ Sp/p and ~ Vp ~ crp ~ e'^. Identifying all 
the terms of order e^™ in equation (Eoh we have 



1 



+ 



+ 



2'"m! 9sj|'" 
1 



kp(s) 



92 



2™-i(m- 1)! asjj'""' 

Q2m-2 

2™-2(m - 2)! dsf:""-^ 



kl>(s)]™-Vp(s) 



[al{s)r-'^r(s) 



;l/|(s) 



(Bl) 



This contains terms involving different powers of s. The terms whose s dependence is s 



2— 7— 2m 



fall off slowest with increasing 



s and these terms dominate at large separations. Retaining only these terms the above expression becomes 



29iF" 



■(^p(s) + 



(B2) 
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Summing up all the terms with different powers of we obtain 
^ a^' d'^ 



9 = II 



where the terms in the square bracket is the linear redshift two-point correlation function (,s(s) (equation | 
up all the terms following a procedure similar to that used in section 2 (equation ^ to and obtain 

Cs(sx,S||) = / G(S||, %/2cr)^f'(sx,S|| +sj|). 



(B3) 

We next sum 
(B4) 
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Figure Bl. This shows the various spherical harmonics of the redshift space two-point correlation function calculated using cq. ( |27 
(NON-LINEAR) and eqs. (||)-(|^) (LINEAR) for a COBE normalised CDM power spectrum with T = 0.25 and h = 0.5. 
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Figure B2. The plots show R2 = £,2/ (,2 different values of /(Q) and 6 for a COBE normalised V = 0.25, h = 0.5 CDM power 
spectrum. 
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